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The proton and neutron structure functions are modified in the nuclear environment [|TJ , 
referred to here as shadowing. In most shadowing models, such as gluon recombination, 



the structure function modifications should be correlated with the local nuclear density. 



The E745 experiment studied the spatial distribution of structure functions with vN 



interactions in emulsion M. They found evidence of a spatial dependence but could 
not determine the form. The spatial dependence of shadowing is reflected in particle 
production as a function of impact parameter, 6, which may be inferred from the total 
transverse energy, Et, produced in a heavy ion collision ||. 

At RHIC and LHC perturbative QCD processes are expected to be an important compo- 
nent of the total particle production. In particular, at early times, Tj ~ 1/pt < 1/Po ~ 0.1 
fm for p ~ 2 GeV, semihard production of low pt minijets sets the stage for further evo- 
lution @. 

The average transverse energy of minijet production is proportional to the initial energy 
density and is the hard scattering contribution to the average total transverse energy. It 
is calculated within a specific detector acceptance e(y). At leading order, the average 
transverse energy for a given minijet flavor, /, as a function of impact parameter is 

E f T (b) = — = Kjct J d 2 rdpTdy 2 dydzdz'^xiFf L (xi,PT,r,z)x 2 Ff(x2 J PT,b - f,z') 
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Identical particles in the final state are accounted for by the factor 1/(1 + Ski)- The ratio 
of the NLO to LO jet cross sections, Kj et , indicates the size of the NLO corrections. A 
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conservative K^ et = 1 || gives a lower limit on minijet production. The cutoff po is the 
lowest pt scale at which perturbative QCD is valid. 

The nuclear parton densities, F t A , are assumed be the product of the nucleon density 
in the nucleus /^(r, z), the nucleon parton density f^(x,Q 2 ), and a shadowing function 
S l (A, x, Q 2 , f, z) where A is the atomic mass number, x is the parton momentum fraction, 
Q 2 is the interaction scale, and f and z are the transverse and longitudinal location of 
the parton so that F^(x,Q 2 ,f,z) = pA(r,z)S l (A,x,Q 2 ,r,z)f^(x,Q 2 ). In the absence 
of nuclear modifications, S l (A, x, Q 2 , r, z) = 1. If the nuclear structure functions are 
homogeneous, then the spatial effects factorize and E T (b, po) is proportional to the number 
of nucleon- nucleon collisions, Tab (b). We have calculated our results with GRV 94 LO 
H and MRST LO parton distributions. The GRV low x gluon density is larger than 
the MRST and produces strikingly different results at LHC while the RHIC results are 
rather similar. We use three different parameterizations of shadowing, all based on nuclear 
deep inelastic scattering. The first parameterization, Si(A,x), treats quarks, antiquarks, 
and gluons identically without Q 2 evolution M. The others evolve with Q 2 and conserve 
baryon number and total momentum. The second, S\{A, x, Q 2 ), modifies the valence 
quarks, sea quarks and gluons separately and includes Q 2 evolution for Q = 2 < Q < 10 
GeV 0. The third, S\(A, x, Q 2 ), evolves each parton type separately above Q = 1.5 GeV 
[ TO . We have studied the spatial dependence of shadowing using two parameterizations. 
If shadowing is proportional to the local density PJLl] , 



S l ws = S\A,x,Q 2 ,f,z) = l + N ws [S\A,x,Q 2 )-l]^^ , (2) 

Po 

where iVws is chosen so that (A/A) J d 2 rdzpA(r, z)S^ s = S\ If instead, shadowing arises 
due to multiple interactions of the incident parton, shadowing should be proportional to 
the path length of the parton through the nucleus, 

SMA.,0 2 ,?,,) = { i + ^(S'(A.,0 2 )-i)Vi-(^) 2 r<R A 

{I r>R A 

The normalization, obtained after averaging over Pa{^i z ), is slightly larger than iV-ws- 
In both cases, at distances much greater than Ra, the nucleons behave as free particles 
while in the center of the nucleus, the modifications are larger than the average value S l . 
We have also studied inhomogeneous shadowing assuming that shadowing is proportional 
to the path length through the nucleus and find that the central shadowing is somewhat 
stronger than with S^ s in Eq. (0). 

The effect of the inhomogeneous shadowing, S^s i n Eq. (0), on E T (b), calculated 
with the GRV 94 LO parton densities, is shown in Fig. [I] for CMS and STAR. The 
CMS and ALICE ratios are similar as are the STAR and PHENIX results. When x 
lies in the shadowing region, central collisions are more shadowed than the average. In 
the antishadowing region, central collisions are more antishadowed than the average. 
When b ~ Ra, the homogeneous and inhomogeneous shadowing are approximately equal. 
When b ~ 2Ra, shadowing/ antishadowing is significantly reduced. As b increases further, 
S — 1 is approached asymptotically. When inhomogeneous shadowing is assumed to be 
proportional to 5*^, Eq. (0), we find stronger shadowing at low impact parameter and a 
somewhat faster approach to S = 1 at large b. 
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Figure 1. The impact parameter dependence of the first Et moment, compared to that 
with S = 1 in CMS, \y max \ < 2.4, at the LHC (left) and in STAR, |y max | < 0.9, at RHIC 
(right) calculated with the GRV 94 LO distributions and po = 2 GeV. The upper plot 
shows the ratio for quarks and antiquarks, the middle plot is the gluon ratio and the lower 
plot is for the total. The horizontal lines show the homogeneous shadowing results: dashed 
for Si, dot-dashed for S 2 , and dotted for S 3 . The inhomogeneous shadowing results using 



S^s f° r Si, circles, S 2 , squares, and S 3 , triangles are also shown [12 



We assume that the transverse energy distribution, da/dE?, can be approximated by 
the Gaussian 

da r d 2 b / [E T -E T {b)} 2 \ 

The total Et distribution is a convolution of the hard and soft components with mean 
and standard deviation 

E T (b) = Y,E T (b)+T AB (b)e al(b)=J2W-^^ + TAB(b)(ei-4) ■ ( 5 ) 

The E T distributions are shown for CMS and STAR in Fig. @. At the LHC, 90% of the 
average Et comes from the hard component, causing the maximum Et to be halved due to 
shadowing. In fact, the low x gluon density is high enough for the incoming gluon to have 
multiple hard collisions, up to 5 in central collisions with no shadowing, as it traverses the 
nucleus when the GRV 94 LO distributions are used. In contrast, an incoming quark has 
at most 1 collision. Similar results with MRST LO distributions typically halve the hard 
Et and significantly reduce gluon multiple scattering. Gluons suffer 1.4-2 hard collisions 



(4) 



4 



while quarks scatter less than once [12] ■ At RHIC, since the hard and soft components 
are comparable, the maximum is shifted by only ~ 7% when shadowing is included. 
Indeed, for 63, the maximum is slightly increased. Thus shadowing significantly affects 
the initial conditions at LHC but not at RHIC. Ref. |]12[ gives more details and discusses 
inhomogeneous shadowing effects on J/ip and Drell-Yan production. 
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Figure 2. The E T distribution predicted for CMS, |y max | < 2.4, at the LHC and STAR, 
ll/maxl < 0-9, at RHIC calculated with the GRV 94 LO distributions and p = 2 GeV. 
The homogeneous shadowing results for all b are: no shadowing (solid), Si (dashed), S2 
(dot-dashed), and S3 (dotted). Modified from Ref. |T2 . 
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